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I. INTRODUCTION 

One of the main goals of microscopic nuclear structure calculations, which are based on realistic models for the 
nucleon-nucleon (NN) interactions, is to obtain a reliable prediction for the equation of state of nuclear matter under 
the extreme conditions, which one has to consider for astrophysical scenarios like a supernova or objects like a neutron 
star. For that purpose one considers NN interactions, which are adjusted to describe the properties of two nucleons 
in the vacuum, i.e. the NN phase shifts, and tries to develop a many-body theory which yields a good description 
for the bulk properties of "normal" nuclear matter, the saturation point of symmetric nuclear matter and properties 
of finite nuclei. A theory, which is able to link such a realistic NN interaction to the properties of nuclear matter at 
normal densities, should provide reliable results for the properties of nuclear matter at higher densities, temperatures 
and isospin-asymmetries, as they occur in the astrophysical objects mentioned above[l|. 

A basic problem of such nuclear structure investigations is the existence of strong tensor short-range components in 
such a realistic NN interaction, which makes it necessary to account for corresponding correlations in the nuclear wave- 
function. In fact, simple Hartree-Fock or mean- field calculations using such realistic NN interactions yield unbound 
nuclei 2]. Rather powerful many-body techniques have been developed to account for correlations beyond the mean- 
field approximation and using rather sophisticated Monte-Carlo techniques one is able to derive the properties of 
light nuclei from such realistic NN interactions Q. In order to obtain a good agreement with the experimental data, 
however, one has to introduce a three-body force. 

An additional three-body force is also required, if one wants to reproduce the saturation point of symmetric nuclear 
matter within a non-relativistic many-body approach based on realistic NN interactions. On the other hand, however, 
it is known already for more than 15 years that the consideration of relativistic features, as it is done e.g. in the Dirac 
Brueckner Hartree Fock (DBHF) approximation yields results for the saturation point of symmetric nuclear matter, 
which are very close to the empirical data[3, |J H, 0|- 

This success is based on the relativistic structure of the NN interaction, a feature which is also present in the 
phenomenology of the Walecka model L 8j. The NN interaction in this model is described in terms of the exchange of 
a scalar meson, a, and a vector meson uj. Calculating the nucleon self-energy, S, from such a meson exchange model 
within a mean- field (Hartree) approximation, one finds that the lo exchange yields a component S°, which transforms 
under a Lorentz transformation like the time-like component of a vector, while the scalar meson exchange yields a 
contribution S s , which transforms like a scalar. Inserting this self energy into the Dirac equation for a nucleon in 
the medium of nuclear matter leads to single-particle energies, which are as small as the empirical value of -50 MeV. 
This small binding effect, however, results from a strong cancellation between the repulsive E° and the attractive E s 
component. As has recently been shown in Q the appearence of these two large and cancelling scalar and vector 
fields in matter is not only a consequence of Dirac phenomenology but has a deeper foundation in the structure of 
the nucleon-nucleon force, i.e. it is intimitely connected to the spin-orbit force in vacuum NN scattering. The scalar 
component S s , which is conveniently described in terms of an effective Dirac mass, leads to a significant enhancement 
of the small component of the Dirac spinors in the nuclear medium. It is this density-dependence of the Dirac spinors, 
which is responsible for the fine-tuning in the nuclear structure calculation, which is necessary to obtain the empirical 
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saturation point. Non-relativistic studies, which cannot account for this feature, may include this effect in terms of a 
three-body force. Moreover, the excitation of anti-nuclcons, automatically included in the relativistic formalism, gives 
rise to a class of three-body forces (Z-graphs) which have to added in the non-relativistic formalism. Doing so, the 
saturation points from non-relativistic Brueckner calculations are shifted towards their relativistic counterparts [fOj . 

A neutron star, however, cannot completely be described in terms of homogeneous nuclear matter in /3-equilibrium 
at various densities. The crust of such a neutron star in particular is a challenge for theoretical nuclear structure 
physics as it incorporates the transition from isolated nuclei via a crystal-like structure of quasi-nuclei embedded in 
a sea of neutrons to phase of homogeneous baryon matter. These intermediate structures have been described by 
means of the Wigner-Seitz (WS) cell approximation employing the Thomas-Fermi approximation based on a purely 
phcnomenological energy-density functional. Such calculations predicted a variety of geometrical structures: Spherical 
quasi-nuclei, which are favored at small densities, merge with increasing density to strings, which then may cluster 
to parallel plates and so on. These geometrical structures have been the origin for the popular name of this phase: 
Pasta phase. 

A step towards a more microscopic study has been made by performing self-consistent Hartree-Fock calculations 
based on effective nucleon-nucleon interactions like the density-dependent Skyrme forces [HI, [12]. Such calculations 
have been performed by Bonche and Vautherin[13j and by a few other groups. These studies show that shell effects 
have a significant influence on details like the proton fraction of the baryonic matter in the inhomogeneous phase [bj]. 
They also provide the basis for a microscopic investigation of pairing phenomena, excitation modes and response 
functions as well as the effects of finite temperature. 

Such self-consistent calculations for inhomogeneous matter have typically been performed assuming a WS cell of 
spherical shape. This geometry, however, does not support the description of triaxial structures, which are typical 
for the pasta phase. Furthermore the transition to homogeneous matter cannot be described in a satisfactory manner 
within such a spherical WS cell [Til . 

Therefore cubic WS cells have recently been considered for Hartree-Fock calculations using effective Hamiltonians 
like the Skyrme forces [3]. These effective Hamiltonians, however, are of purely phenomenological origin. They are 
adjusted to describe nuclear matter at normal density and finite nuclei. Therefore one cannot expect any predictive 
power of such calculations if one extends its region of application to higher densities, large proton-neutron asymmetries 
and finite temperatures. 

Therefore in the present work we will perform relativistic mean field calculations employing a model for the La- 
grangian, which is based on microscopic DBHF calculations using realistic NN interactions. Attempts to derive an 
effective Lagrange density by fitting density-dependent coupling constants in such a way that a mean-field calculation 
reproduces details of a DBHF calculation have been made before[l(|[l7|- The calculations presented here are based on 
DBHF calculations of van Dalen et al.[l|j| for asymmetric nuclear matter. These results for the DBHF calculations for 
homogeneous asymmetric nuclear matter have successfully been used in bulk studies of neutron stars [19(. Therefore 
we have tried to ensure that the effective mean field parameterization reproduces the properties of the equation of 
state derived from this DBHF calculations at high densities with good accuracy. Furthermore we introduce a small 
correction term, which ensures that the mean-field calculations yield a good description of binding energy and radii 
of finite nuclei as well. In this way we obtain an effective field theory, which is based on a realistic model for the NN 
interaction in the vacuum, reproduces the details of microscopic DBHF calculations for asymmetric nuclear matter 
at high densities and yields good agreement with the empirical data for the saturation point of asymmetric nuclear 
matter as well as bulk properties of finite nuclei. Therefore the resulting parameterization should be a good candidate 
to determine a reliable equation of state for baryonic matter in a large interval of densities. 

Using this Lagrangian we will present results of relativistic mean-field calculations in Cartesian WS cells considering 
a range of densities in which the transition from isolated nuclei to homogeneous matter occurs. Special attention will 
be paid to the effects of finite temperature on the formation of the geometrical structures. 

After this introduction section 2 contains a short review of the Density Dependent Relativistic Mean Field (DDRMF) 
approach and its application to the description of infinite matter as well as finite nuclear systems and baryonic 
structures in a Cartesian WS cell. The parameterization of the density dependent coupling constants and the fit to 
the DBHF results is described in section 3. Results for the inhomogeneous structures, which are typical for the crust 
of neutron stars are presented in section 3. Special emphasis is made to explore the effects of finite temperature and 
the differences as compared to studies using a purely phenomenological Skyrme forces. 

II. THE DDRMF APPROACH 

The Density Dependent Relativistic Mean Field (DDRMF) approach is an effective Field theory of interacting 
mesons and nucleons. Following the usual notation we consider scalar (<j, 5) and vector mesons (u),p), which with 
respect to the isospin correspond to isoscalar (<r, to) and isovector (5, p), respectively. The Lagrangian density consists 
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of three parts: the free baryon Lagrangian density £b, the free meson Lagrangian density Cm and the interaction 
Lagrangian density C mt : 



C — Cb + Cm + Ant) (1) 



which take the explicit form 



-5 E (f f W^^r m ^W^( K )j. ( 2 ) 

K—uj.p.'y 

Ant = - 5<T** ff * - .95*t-* 5 * 

with the field strength tensor -F( K w = 9fj,Ar K \ v — d v A^ K ^ for the vector mesons. In the above Lagrangian density the 
nucleon field consisting of Dirac-spinors in isospin space is denoted by \P and the nucleon rest mass by M — 938.9 
MeV. The scalar meson fields are $<j and $>g, the vector meson fields A^) and Ar p y Bold symbols denote vectors 
in the isospin space acting between the two species of nucleons. The mesons have rest masses m K for each meson k 
and couple to the nucleons with the strength of the coupling constants g K . The electromagnetic field couples to 
the nucleons by the electron charge e 2 = Ana where a is the fine structure constant. Notations are taken from [20| : 
x = and x p denote the contravariant and covariant vectors in space-time, 7 M the Dirac 7 matrices and r consists 
of the isospin Pauli matrices r^. 

Allowing for a density dependence of the baryon-meson vertices improves the capability of the model significantly, 
since it enables the latter to assimilate the self-energies of a Dirac-Brueckner-Hartree-Fock calculation. This means 
that the coupling constants g K depend on a density p(*ff, ^) obtained from the nucleon field 'J. In the literature we 
find a dependence on the scalar density or on the vector density [21] . It turned out that the dependence on the zero 
component of the vector density, the baryon density p = Tr(^7o^ / ), is the most suitable one since it describes finite 
nuclei better and has a natural connection to the vertices in the DBHF calculations [l7| ■ The functions for the various 
couplings are specified later together with the associated parameter set. 

Applying the variational principle to the Lagrangian we obtain a Dirac equation for the nucleons and Klein-Gordon 
and Proca equations for the meson fields [22| • 

Due to density dependent vertices the variation principle changes to 

SC = dC B£5p 

3* dp <S§ ' U 

where the second expression creates the so called rearrangement contribution E# to the self-energies of the nucleon 
field. These rearrangement contributions contribute only to the zero component of the vector self-energy. Including 
these additional contributions we denote the Dirac equation for the nucleonic single-particle wave function ip a in 
Hartree approximation 

(up + (So + E fl ) + (3(M + E<?)) tp a = e a (4) 
where the self-energy contributions read 

E s = g a $ a +55$5T 3 , 

So - g u A^ + g p A { P) T3 + e l -{\ - r 3 )4 7) , (5) 
and the rearrangement self-energy contribution E# is obtained by 

v f d 9°<f, s , d 9s , s . 9g u (u) ) dg p (p) \ 
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The various densities are obtained from the nucleon single-particle wave functions in the "no-sea" approximation 

as 



p s {x) 


Oi 




Oi 


p{x) 


= y^?7a ^q(x)7Q'0 a (x) 
a 


Pa(x) 


= ^ T) a jl a (x)^ T 3 ^ a (x) 

Oi 


p {em) (x) 


= y^T?a^q(cc)|(l - T 3 )lp a (x) [-p e (x 

Oi 



(7) 

where p s is the scalar density, p the baryon density, p% the scalar isovector density, pz the vector isovector density, and 
p{em) charge density. The occupation factors rj a have to be determined from the desired scheme of occupation. 

Neglecting retardation effects the Klein-Gordon equations reduce to inhomogeneous Hclmholtz equations with 
source terms 22] 

(-A + ml) $ CT = -g„ p s 

(-A + TT#)$ 4 = -95 Ps 

(-A + m»)4 w) = 9»P 
(-A + m 2 p )A { p) = g pP3 

-A4 7) = ep [em \ (8) 

from which the self-energy contributions ([5|) are obtained. The Dirac equation for the nucleons (jlj, the evalua- 
tion of the resulting densities 0, these meson field equations © and the calculation of the resulting self-energy 
contributions ([5]) form a set of equations, which have to be solved in a self-consistent way. 



A. Nuclear Matter 



In nuclear matter the electromagnetic field is neglected and translational invariance is assumed. The Dirac spinors 
u(k, s, i) with momentum k, spin s and isospin i are solutions of the following Dirac equation: 

[7 • k* + m*]u(k, s, i) — j n E*u(k, s, i), (9) 

where the starred quantities are defined by 

to* = M + 

E* = £(fc)-£ ,i(fc), (10) 

and the on-shell condition reads 

Ef=k 2 + mf. (11) 

The positive energy solutions of the Dirac equation ^ are obtained in straight analogy to the results for the free case 
by 

E* + to* \ 1/2 / 1 



u{k, s,i)=[ \ 2Ef 1 ) I cr.kT ) xi/a(g)xi/2(Q, (12) 



BT- 



where the spinors Xi/2(s) and Xi/2(i) account for spin and isospin projection. 

From these spinors the vector and scalar densities for protons and neutrons can easily be calculated in infinite 
matter as 



j m{k)k 2 dk 



Ei(k)* 



Pi = \ I m(k)k 2 dk. (13) 
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The energy density £ and the pressure P are obtained from the energy-momentum tensor 



t—p,n 



+lj2(^s,iPt + ^o,iPi) (14) 

p=lh Tjj ) = fc4 



j=i i=p,n yk 2 + m *C 

+Z RP +lj2 p\ + £o,< Pi) • (15) 

Rearrangement does not contribute to the energy density but it modifies the pressure, which coincides with the 
corresponding thermodynamical relation. 
The entropy density S is obtained as 

S = E / ^( fc ) fc2 dk ^( fc ) ln ^( fc )) + C 1 - /*(*)) ln (! - Mk))} , (16) 

i—p,n 

where fi(k) are the thermal occupation probabilities. Note that Tji and fi may differ, e.g. if pairing is introduced in 
the finite temperature BCS formalism. The energy density £ and the entropy density S yield the free energy density 

T = £-TS (17) 

and the thermodynamic potential 

= £ - TS - p p Pp - fJ, n Pn- (18) 

Minimizing the thermodynamic potential Q for fixed meson fields we obtain the chemical potentials 

Pi = £f,i + Eo,i + £r (19) 
from the Fermi energy ep,i = E*(k F ,i) an d the thermal occupation factors in the "no-sea" approximation 

m = i + cM(E;-e F , t) /Ty i = p ^ (20) 

The Fermi momenta k F ,i for protons and neutrons have to be determined from the given proton and neutron densities 
(see eq. [T3l) . A self-consistent calculation including densities, Fermi momenta, effective masses and the self-energy 
contributions for protons and neutrons yields all described quantities. 

B. Finite Systems 

For finite nuclei and the description of nuclear matter in a Wigner-Seitz cell (WS) the Dirac equation (0| and the 
meson field equations ([8]) are solved in spatial representation. The numerical procedure to solve the Dirac equation in 
the cubic box is the same as in [l5j with an extension of the code to allow for the density-dependent coupling constants. 
Pairing correlations are included like in flil . [l5j . To be able to perform consistent finite temperature calculations the 
finite temperature BCS formalism (FT-BCS) has been included in the code which allows the simultaneous evaluation 
of pairing and finite temperature effects. 

According to Goodman [23[ the normal and anomalous occupation factors of the Hartree-Fock single particle states 
are obtained in the FT-BCS approach as 

Va = (l-2f a )v 2 a + f a , 

Ca = (l-2f a )u a V a , (21) 

where v a and u a denote the usual quasi-particle occupation factors and f a the thermal occupation factors 

fa = 1 + exp^/T] (22) 
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depending on the quasi-particle energy E a — yj (e a — £F,i) 2 + A^ , where e a are the single-particle energies and A a 
the state-depending pairing gap. 

The normal occupation r\ modifies the densities (|7|) and (|13j) , while the anomalous occupation factors enter into the 
anomalous density 



(r) = i^CakWr)| 2 . (23) 



X 

The local gap function based on a density dependent pairing force V of zero range 

A(r) = -V» X (r), (24) 

modifies the state-dependent pairing gaps 

A Q = J d 3 r A(r) |^ Q (r)| 2 . (25) 

These pairing gaps A a are evaluated in a self-consistent procedure fixing the Fermi energies for protons and neutrons 
Ep t i by the corresponding particle number conditions 

# = 5>a. (26) 
a 

Finally the the pairing energy E pa i r is obtained from the state dependent pairing gaps 

Epair — ^ ^ ] A a (^ Q . (27) 

a 

The center of mass correction which gives a significant contribution to the binding energy of light nuclei is treated 
in the usual harmonic oscillator approximation like in [17| 

E cm = ~hu (28) 

with Huj = 41 A~ 1 / 3 MeV. 

The total ground state energy includes the center of mass correction, the pairing energy and the energy of the 
relativistic mean field 

Eq = EftMF + Ep air + E cm , (29) 

where the Hartree ground state energy is obtained from the Dirac Hamiltonian corresponding to the Lagrangian 
density (JlJ similar to [22j . Considering the rearrangement and the single-particle energies the ground state energy is 
calculated by 



Ermf = ^2 r l a£a ~ I d3r s «( r ) P( r ) 

a J 

~\ E Jd 3 r[Es,i(r)pKr) + ^oAr)Pi(r)} 



(30) 



III. DENSITY DEPENDENT PARAMETERIZATION FROM DBHF THEORY 



In this section we will describe the parameterization of density dependent coupling constants to be used in the 
DDRMF approach, which is based on the DBHF calculations for homogeneous asymmetric nuclear matter of van 
Dalen et al.fla]. which employ the meson exchange interaction model OBEPA of the Bonn potentials The aim is 
to obtain a parameterization, which accurately reproduces the results of the DBHF calculations at high densities but 
also provides a good description of bulk properties of finite nuclei. This should lead to an equation of state covering 
a very broad range of densities as e.g. in Refs. [3, HE]- 

A first attempt to translate the results of DBHF calculations in terms of density-dependent coupling constants 
would be to consider the scalar (S s ,i) an d vector contributions (Sq^) for protons and neutrons calculated in the 
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DBHF approximation for asymmetric nuclear matter of a density p and proton - neutron asymmetry p$ and equate 
these self-energy terms with the corresponding mean-field expressions. This leads to the following expressions for the 
effective coupling constants for the a, lj, S, and p mesons [35| : 



9cr(p, P3) 

9uj(P,P3) 

9s(p,P3) 
ms 





t ^s,n \^Fn ) 


O 
Z 


P 


1 ^0,p(^Fp) 


+ £o, n (fcFri) 


2 




1 S S! j,(fc Fp ) 


— S Sin (fcF„) 


2 


P3 


1 S 0:P (fcFp) 


— Xo.ri^Fn) 


2 


P3 



(31) 
(32) 

(33) 



( 9p{P,P3) \ 2 1 So,p(fc F p) ~ Eo,„(fcF«) , . 



with p s = p s p + p s n , p = p p + p n , p s 3 = p s p - p s n , and p 3 = p p - p n , where 



1 = (35) 



^ ^0 ' y/mf + k 



and 



d 3 fc - H (36) 



(2tt)3 y 3^ 

are, respectively, the scalar and vector density of particle (i = n,p). In contrast to widely used RMF theories we 
explicitly include the scalar isovector meson 5 since this provides a mechanism to account for the differences in the 
scalar self-ener gies and the corresponding effective Dirac masses for protons and neutrons in asymmetric matter [10]. 
Following Ref. [26[ the inclusion of the S meson in the DDRMF theory has important consequences for the dynamics 
of neutron-rich nuclei. In addition, this meson is important for astrophysical applications since the dense asymmetric 
matter gets softer, which means that the pressure rises more slowly for larger densities. 

The expressions of eqs. (|3ip - (|M|) identify coupling constants, which depend on the density p and the proton - 
neutron asymmetry (^3) Since, however, the dependence of the coupling functions turns out to be weak (see e.g. 
[3[lM[13)> we have ignored this dependence on p 3 to keep the DDRMF functional as simple as possible. 

The fact that a kind of renormalization is required when DBHF results are mapped on Relativistic Mean Field 
(RMF) theory has already been pointed out in Ref. [l?], Gil • The reason is that two essential differences exist between 
DBHF and RMF concerning the structure of the self energy. First, the DBHF self energy terms explicitly depend on 
the momentum of the particle, a feature which is absent in RMF. This reflects the non-locality of the DBHF self-energy 
terms, which originates from the Fock exchange terms but also from non-localities in the underlying NN interaction. 
Since this explicit momentum dependence of the DBHF self energy components is generally weak [2?], [28| , we will 
ignore this non-locality effects and always consider the DBHF self-energy terms calculated at the corresponding Fermi 
momenta kpi as already expressed in eqs. (|3"T1) - (|M1) . 

Secondly, the appearance of a spatial contribution of the vector self energy XV in the DBHF theory, which is not 
present in the RMF model. The Sy component originates from Fock exchange contributions which are not present 
in the RMF theory. For an accurate reproduction of the DBHF energy functional the spatial Ey component has to 
be included in a proper way. The effects of the the Ey component in the Dirac equation for homogeneous nuclear 
matter can be absorbed into a renormalization of the scalar and time like vector component of the self-energy. This 
leads to an effective Dirac mass, which has to be identified with the RMF effective mass, i.e. 

_ M + T, S!i (k Fi ) _ RMF 



1 1 + 

This leads to the renormalized scalar self energy component 



AZ + Eff. (37) 



^RMF _ £s,i( fcgi) ~ ME Vti ( k Fi ) 
1 + E v ,i(fcfi) 



•^RMF _ ^s,i\^t%} ^^v : t\^^i/ 



In a corresponding way the following expression for the renormalized vector self energy component is obtained 

^rmf v (j. \ Z v ,i(k Fi )[3E F iPi + m*pf\ 

h 0,i =^0,U«fiJ 1— ■ (39) 
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T P T 
J 1 


m [MeV J a K o K c K d K 


a 


0+ 


550 7.7868 2.58637 2.32431 3.11504 




1" 


782.6 9.73684 2.26377 7.05897 


8 


0+ 1 


983 0.503759 2.68849 6.7193 0.403927 


P 


1" 1 


769 4.56919 5.45085 1.20926 



TABLE I: Parameter set from DBHF by van Dalen et al. [3 for the density dependent relativistic mean field approach. 





DBHF DDRMF 


Paat [MeV] 
E/A [MeV] 
K [MeV] 
a s [MeV] 


0.181 0.166 
-16.15 -16.23 
230 335 
34.36 29.92 



TABLE II: Nuclear matter properties obtained in the different models. The DDRMF model includes the correction from 
renormalization. 



These renormalized self-energy components are now inserted into eqs. ([3lTl -(l34 p to obtain the renormalized density 
dependent coupling functions. 

The density dependent isoscalar couplings are obtained from eqs. (|3l"1) - (j34|) using the symmetric nuclear data 
((3 — 0.0) up to a density of 0.5 fm -3 , whereas the isovector ones are obtained using the neutron matter data of the 
DBHF calculations in Ref. [ufy. 

In order to make this parameterization easily accessible, we have fixed the masses of the mesons to the corresponding 
masses in the underlying free NN interaction (see table [TJ> and parameterized the density dependence of the coupling 
constants by 



[b K + d K x 3 ] exp(-c K a;), 



(40) 



where x = ps/poi an d po = 0.16 fm -3 is the saturation density of symmetric nuclear matter. The values obtained for 
the fit of the coupling functions are summarized in table |TJ 

Employing these functions in a DDRMF calculation of finite nuclei, we obtained binding energies and radii, which 
are too small as compared to the empirical data. This is in line with the observation that also the underlying DBHF 
calculations for symmetric nuclear matter yield a saturation density, which is too large as compared to the empirical 
result (see table [TTJ) . In order to cure these problems, we have considered a slight reduction of the coupling constant 
for the u> meson at densities around the saturation density in a form: 



•{Pb) = 9u>{p) - a cor exp ( - 



Po 



(41) 



Adjusting the parameters of this correction to a cor = 0.014 and b cor = 0.035 fm -3 , we obtain an improvement for the 
saturation density of nuclear matter (see table [TTJ) a rather good description for energies and radii of finite nuclei (see 
table ITTT)) . The lighter nuclei are a little too much bound but the binding energies for heavy nuclei are well reproduced. 
The charge radii r c show a good agreement for the lighter nuclei whereas for the heavy nuclei the radii are a little too 
small. 





lb O m Ca m Ca ' M Zr Pb 


E/A [MeV] 
E/A exp. [MeV] 


— 8.35 -8.73 -8.73 -8.74 -7.87 
-7.98 -8.55 -8.67 -8.71 -7.87 


r c [fm] 
r c exp. [fm] 


2.78 3.44 3.45 4.17 5.31 
2.74 3.48 3.47 4.27 5.50 



TABLE III: Results for closed shell nuclei applying the DDRMF parameterization from the DBHF results. Experimental values 
are taken from [l~7l ] . 

In table |TT] the nuclear matter properties of the DBHF calculations are compared to those of the DDRMF using 
the adjustment of g u according to ea. (|41[) . Due to this adjustment the nuclear matter properties are slightly changed 
around the saturation point. The saturation density is shifted to a lower value so it is even closer to the experimental 
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value. The energy per nucleon E/A of symmetric nuclear matter is fairly well reproduced but the compression 
modulus K is rising due to the fit procedure. In this context it is worth to mention that RMF fits to finite nuclei 
require relatively high compression modulus K of about 300 MeV [29| ■ EoSs with a stiff high density behavior stand, 
however, in contrast to the information extracted from heavy ion reactions 30]. Note, however, that the correction 
for the u> coupling function is restricted to densities around po and should not affect the high density behavior of the 
EoS. 

The symmetry energy coefficient a s is decreased by about 4 MeV compared to the DBHF results. This reduction of 
the symmetry energy is connected to the decrease of the saturation density as the two values are read of at essentially 
different densities. These changes should not affect the astrophysical properties of the model in the high density 
region, what is important to be able to reproduce the heavy neutron stars. Moreover, the value of a s ~ 30 MeV at 
Psat — 0.166 fm~ 3 agrees well with the presently favored experimental value of a s ~ 31 ± 1 MeV (3l| . 

IV. RESULTS AND DISCUSSION 

In this section we are going to discuss results of density dependent relativistic mean field calculations (DDRMF) 
employing the parameterization of effective coupling constants as described in subsection . Pairing correlations 
are included in terms of the BCS approximation assuming a density dependent zero-range pairing force, which has 
already been used in earlier investigations [TJ, [l5[ • 

The calculations are performed in a cubic Wigner Seitz (WS) cell. The origin of the coordinate system is put in 
the center of the box and we assume the density profiles to be symmetric under reflection on the x = 0, y = and 
z = planes. Therefore we can restrict the calculation to a cubic box of length R in each direction, so that the "size" 
of this box R in each direction is half the length of the WS cell. This box size R has been adjusted to minimize the 
total energy per nucleon for the density under consideration. 

All calculations which we discuss in this section for charge neutral matter containing protons, electrons and neutrons 
in /3-equilibrium. This implies that the chemical potentials of these species obey 

A*n = Mp + Me- (42) 

We start our discussion by exploring the existence of geometrical structures in the density distribution of protons 
and neutrons, which are typical for the "pasta-phase" in the crust of neutron stars. For that purpose we display such 
density distributions obtained from DDRMF calculations at zero temperature in the panels on the left hand side of 
Fig. [T] 

The top panel in the left column represents a nuclear structure at a baryonic density of 0.020 fm 3 at zero temper- 
ature. In this case the density profiles are identical in all three Cartesian directions. In the center a quasi-nuclear of 
the WS cell droplet structure is formed with a neutron sea between these quasi-nuclei. 

Increasing the baryonic density to p = 0.057 fm~ 3 , we obtain the second panel in the left column of Fig. Q] In this 
case the proton density shows almost the same value along one axis, which is chosen to be the z-direction (dashed 
curves), while it drops to zero, if the values for x and y tend towards R. Here we find a so-called rod or Spaghetti- 
shape structure for the proton density distribution. Note, however, that the density distribution is not as simple as 
these names rod-structure or Spaghetti-structure suggest. This can be seen from the complete density-distribution in 
the y = (or x — 0, which is identical for this geometry), as it is displayed in Fig. [21 There is a dip in the center 
of the quasi-nuclear structure and the density distributions exhibit a structure, which is close to those of nuclei with 
strong prolate deformations, which are aligned and touch each other along the z-axis. 

At larger densities the density distributions show again the same distribution along the axis in all Cartesian 
directions. It is still slightly enhanced in the center and drops only by a small amount from the center to the 
boundary of the cell. An example of such a structure at the baryon density of 0.067 fm~ 3 is displayed in Fig. [T]in the 
bottom panel of the left column. However, the proton density distribution in the z = plane displayed on the left in 
Fig. [3] provides a more detailed picture of the geometric structures at this density: The proton density drops along the 
diagonals, which implies that the regions of high densities in the center of the WS cells are connected by arms along 
the Cartesian axes. From this point of view this structure may be called a grid structure. If, however, the coordinates 
are shifted in such a way, that the centers of the quasi-nuclei are in the corners of the box as it is displayed in Fig. [3] 
on the right, one finds that this structure can also be interpreted as a bubble structure. Summarizing we observe a 
grid structure with bubbles in between, which provides a last step of a smooth transition to homogeneous matter. 

These kind of geometrical structures and others have also been obtained in a recent investigation employing the 
Skyrme Hartree-Fock approach [l5| . Some typical results for the density profiles at zero temperature are presented in 
the left column of Fig. Q] These examples represent a droplet structure (p = 0.027 fm -3 , top panel), a rod-structure 
with a reduced density at z = R as compared to z = {p — 0.0565 fm -3 , middle panel) and a rod-structure with 
constant densities along the z-axis (p = 0.0625 fm -3 , bottom panel) 
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FIG. 1: (Color online) Density distributions resulting from density dependent relativistic mean field (DDRMF) 
calculations for protons (black color) and neutrons (red color) as a function of Cartesian coordinates x, y, z. The 
panels in the left column refer to zero temperature calculations at the densities 0.020 fm -3 (top), 0.057 fm~ 3 , and 
0.067 fm~ 3 (bottom), while those in the right column are obtained for the temperature T — 5 MeV at the densities 

0.018 fin" 3 (top), 0.044 far 3 , and 0.057 fm" 3 . 




FIG. 2: Proton density distribution in the z = and the y — plane for the density dependent relativistic mean 
field (DDRMF) calculation at an average density of p = 0.057 fm~3. 




FIG. 3: Proton density distribution in the z = plane for the density dependent relativistic mean field (DDRMF) 
calculation at an average density of p = 0.067 fm~3. The surface plot on the left shows the density distribution with 
the origin in the center, while in one on the right hand side the origin is shifted to the corner. 
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FIG. 4: (Color online) Comparison of density distributions at finite temperature with those at zero temperature 
resulting from Skyrme HF calculations for protons (black color) and neutrons (red color) as a function of Cartesian 

coordinates x,y,z. The panels in the left column refer to T — MeV calculations at the average densities 0.0271 
fm~ 3 (top), 0.0565 fm -3 , and 0.0625 fm -3 (bottom), while those in the right column are obtained for T = 5 MeV at 
the baryon densities 0.0273 fm" 3 (top), 0.0545 fm~ 3 , and 0.0648 fm" 3 . 





Skyrme 
HF TF 


DDRMF 
H TF 


Skyrme 
HF TF 


DDRMF 
H TF 


droplet-rod 
rod-slab 
slab-homogeneous 


0.042 0.066 
0.070 0.078 
0.080 0.085 


0.057 0.048 
0.064 0.061 


0.040 0.048 
0.065 0.048 


0.044 0.045 
0.044 0.045 




T = 


T = 5 MeV 



TABLE IV: Comparison of densities at which shape transitions occur using the Skyrme and the density dependent relativistic 
mean field (DDRMF) approach. Results are compared, employing the microscopic Hartree-Fock (HF), the mean-field or 
Hartree (H) and the Thomas-Fermi (TF) approximation. Results displayed in the left part of the table refer to a temperature 
of T = 0, whereas those listed in the right part were obtained for T — 5 MeV. All entries are presented in fm~ 3 . 



Comparing the DDRMF results of Fig. [T] with these density distributions obtained from non-relativistic Skyrme 
calculations, we see that in case of Skyrme Hartree-Fock calculations the different structures are more pronounced. 
Also we observe a slab structure in the Skyrme Hartree-Fock approach, which is absent in the DDRMF calculations. 
This differences may originate from the finite range of the NN interaction in the DDRMF approach as compared to 
the zero-range forces of the Skyrme model. 

In both cases we find that the microscopic calculations in a Cartesian WS cell for densities in the range of 0.01 fm -3 
to 0.1 fm~ 3 lead to quite a variety of shapes and quasi-nuclear structures with smooth transitions in between, which 
may be characterized as quasi-nuclei, rod structures, slab structures, which are all embedded in a sea of neutrons and, 
finally, the homogeneous matter. In the left part of table HVl the densities, at which the transitions from one shape to 
the other occur at T = 0, are listed for Skyrme HF calculations as well as for the density dependent relativistic mean 
field approach. 

Performing DDRMF calculations at the finite temperature T = 5 MeV we obtain for neutrons and protons density 
distributions as displayed in the right column in Fig. [TJ Inspecting these plots, it is clear that the geometrical 
structures observed for T = persist to some extent also at finite temperatures as large as T — 5 MeV. It is also 
obvious, however, that finite temperature effects tend to dissolve these structures: The structures in the right column 
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FIG. 5: Proton-abundances and energy per nucleon resulting from density dependent relativistic mean-field 
calculations (DDRMF) at different densities. The results evaluated in cubic Wigner Seitz cells (circles) are compared 
to those of homogeneous infinite matter (solid lines) and corresponding Thomas-Fermi calculations (dashed lines). 

at corresponding or even lower densities are less pronounced than those in the left column of Fig. [1] The same features 
can also be observed for the Skyrme Hartree-Fock approach, as can be seen from Fig. @] 

The transition densities for all calculations in the WS cell at a temperature of T = 5 MeV are summarized in the 
right part table IIVI If we compare the Skyrme HF transition densities at T — 5 MeV with the zero temperature 
results, the transition densities are lower at finite temperature and the slab structure disappears. The results obtained 
from the DDRMF approach do not any more show a rod structure, but the transition from the droplet structure to 
homogeneous matter occurs via a grid structure, which can be regarded as transition to the homogeneous phase like 
in the zero temperature case. 

A comparison of the proton abundances and the energy per nucleon obtained within the DDRMF model in a 
Wigner-Seitz cell for zero temperature are displayed in Fig. [5j The solid lines represents the homogeneous matter 
results, while the circles display the results obtained from the microscopic density dependent relativistic mean field 
calculation in the cubic box. At densities larger than about 0.08 fin -3 the proton abundances and energies coincide 
for homogeneous matter and the microscopic calculation, which implies that for densities larger than 0.08 fm -3 the 
formation of inhomogeneous structures does not provide any gain in energy for baryonic matter in /3-equilibrium. For 
lower densities, however, the non-homogeneous structures provide a gain in binding energy up to 1.7 MeV per nucleon. 
These quasi-nuclear structures, embedded in a sea of neutron matter, are also responsible for the enhancement of 
proton abundances at small global densities: the protons are constrained to the regions of enhanced densities. 

One can try to simulate these features resulting from the non-homogeneous structures using the Thomas-Fermi 
approximation. Like in a Ref. [l5| we want to compare the DDRMF results with the Thomas-Fermi calculation based 
on a local density approximation of the energy density for homogeneous nuclear matter. In our Thomas-Fermi (TF) 
calculations we use a parameterization for the density distributions like in 32] 



where the central density pf 1 , the peripheral density p° ut , the structure radius Ri and an exponent U are the variational 
parameters completed by the box-size R and the proton abundance. As an alternative a Wood-Saxon type density 
parameterization has been considered, but it turned out that the results are comparable and only the surface energy 




(43) 
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FIG. 6: (Color online) Comparison of proton abundances and energy per nucleon obtained from DDRMF (black 

color) and Skyrme Hartree-Fock calculations (red color). Solid lines are obtained from a polynomial fit to the 
DDRMF or the Skyrme Hartree-Fock calculation in the cubic Wigner-Seitz cell, while dashed lines represent the 

results for homogeneous infinite nuclear matter. 



constant Fq has to be readjusted. For the description of rod-shape quasi-nuclear structures cylindrical coordinates 
are used to parameterize the dependence of the densities on the radial coordinate in a way corresponding to eq. (|43[) . In 
the case of quasi-nuclear structures in form of slab-shapes these parameterizations are considered for the dependence 
of the densities on the Cartesian z-coordinate. 

Assuming those density distributions, the TF energy is calculated as a sum of the bulk energy, i.e. the integrated 
nuclear-matter energy densities, the Coulomb energy, plus the contribution of a surface term of the form [32l |33| 

E smi = F f d 3 r | Vp\ 2 . (44) 

./WS-ccll 

The parameters of the density distribution in (|43[) are varied to minimize the energy of the system under consideration. 
The Parameter Fq for the surface energy term in (|44| has been adjusted in such a way that the TF calculation 
reproduces the experimental binding energy of the nucleus 208 Pb, which has also been reproduced in the DDRMF 
calculation. This adjustment leads to a value for Fo of 61.0 MeV fm 5 , which is rather similar to the value of 68.3 
MeV fm 5 , derived in [l5[ for Skyrme HF case. 

Such TF calculations lead to non-homogeneous structures similar to those calculated in DDRMF and we obtain 
transitions between the various shapes at the densities, which are listed in table IIVI Comparing the Thomas-Fermi 
results to the one obtained from the microscopic DDRMF calculations in Fig. [5] in the density range below 0.05 fm -3 
it is observed that the proton abundance from the Thomas-Fermi calculations are lower and the energy is larger as 
obtained from the corresponding microscopic DDRMF results. The surface energy constant Fq may be readjusted to 
reproduce the microscopic DDRMF results for the nuclear structures in /3-equilibrium. This leads to a value of Fq, 
which is only about half the value for Fo derived from the fit to the properties of 208 Pb. This suggests to consider a 
surface constant Fq, which depends on the isospin asymmetry of the system. This feature is very similar to the one 
observed applying the Skyrme Hartree-Fock approach [l5T|. 

To enable further explorations on the differences between the predictions for the "pasta-phase" derived from 
DDRMF and Skyrme HF calculations, results of those two approaches for the energy per nucleon and the pro- 
ton abundances are directly compared in Fig [6] The results of both models show a similar overall behavior and both 
predict a transition to the homogeneous phase at a density of about 0.08 fm -3 . Looking into details, however, one finds 
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FIG. 7: Proton abundances and energy per nucleon resulting from density dependent relativistic mean field 
calculations at the temperature T = 5 MeV. The results evaluated in cubic Wigner-Seitz cells (circles) are compared 

to those of homogeneous infinite matter (solid lines). 

rather interesting differences: The energy per nucleon is larger for the Skyrme calculation at all densities considered 
in Fig[6j This is in line with the fact that the Skyrme force SLy4, which we have used here, yields a parameter for the 
symmetry energy a$ of 32 MeV[34j|, which is larger than the one derived from DDRMF (see table [TTJ) . Note, however, 
that the symmetry energy coefficient is determined at saturation density, whereas here we are comparing structures 
at half the saturation density and below. A larger symmetry energy at those small density would also explain the 
larger proton abundances obtained in the Skyrme approach. 

The energy gain due to the formation of non-homogeneous structures is slightly larger for the Skyrme approach at 
medium densities of 0.02 to 0.08 fm -3 . This is in line with our observation made in the discussion above, that the 
geometrical structures determined in Skyrme HF are typically more pronounced than those extracted from DDRMF. 
At very small densities the DDRMF leads to a larger gain in energy forming non-homogeneous matter distributions. 
At these densities the DDRMF calculations in the WS cell also predict a larger proton abundance than the Skyrme 
approach. 

The parameters of the Skyrme interaction have only been adjusted to the data of finite nuclei and the saturation 
point of symmetric matter, i.e. data of nuclear systems at normal densities and proton-neutron asymmetries. The 
parameterization of the DDRMF on the other side also reproduces those data but is furthermore based on a realistic 
NN interaction. Therefore, we trust in the predictions of the DDRMF approach for this region of small densities and 
large proton-neutron abundances more than in those based on the Skyrme Hamiltonian. 

Finally, we are going to discuss the effects of finite temperature. For that purpose Fig. [7J displays results for proton 
abundances and the energy per nucleon evaluated in the DDRMF approach at a temperature of T — 5 MeV. The lower 
panel in Fig. [7J shows two curves: the energy per nucleon (upper curve) and the free energy per nucleon (lower curve). 
Comparing with the results obtained from zero temperature DDRMF calculations we recognize that the results for 
T = 5 show a minimum for the energy per nucleon at a density of about 0.05 fm -3 , while the zero temperature results 
show no minimum. However, the free energy, which is the energy of consideration, is lower than the energy in the 
T = case and rises with a larger slope, which means that the finite temperature increases the pressure derived from 
the free energy. 

The energies resulting from the DDRMF calculations, which allow for non-homogeneous structure yield a gain in 
energy for densities below 0.05 fm -3 . This means that a temperature of 5 MeV lowers the critical density for the 
formation of local structures from 0.08 fm~ 3 at T = to 0.05 fm~ 3 at T = 5 MeV. At values below this critical density, 
the gain in the free energy at finite temperature turns out to be considerably smaller than in the zero temperature 
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limit. This reflects the feature that finite temperature effects tend to suppress the formation of density fluctuations. 
Note, that these small differences between the free energy of the homogeneous and inhomogeneous structures is to 
quite some extent due to the enhancement of the entropy in the inhomogeneous case. We observe the same features 
in Skyrme HF calculations at finite temperature. Finite temperature also reduces the variety of geometrical shapes 
in the Skyrme HF as well as in the DDRMF approach (see table IIV[) . 

Fig. [7] also displays results from Thomas-Fermi calculations. In this case the temperature dependence is contained 
only in the the bulk contribution, i.e. the integrated nuclear-matter energy densities calculated at the finite temper- 
ature. The parameter for the surface energy Fq is still adjusted to reproduce the properties of 208 Pb at T — 0. As 
discussed above one could now try to fit this constant Fq in such a way that we reproduce results for the energy, 
the free energy and proton abundances of the DDRMF approach. This readjustment leads to a value of 27 MeV fm 5 
for the parameterization of ea. (|15)) . The value is lower by about 15% compared to the corresponding one at zero 
temperature. Hence a slight reduction of the surface energy constant Fq could improve the Thomas-Fermi approach 
at finite temperature. 

V. CONCLUSION 

Density dependent relativistic mean field (DDRMF) calculations have been performed to study the structure of 
baryonic matter in /3-equilibrium with electrons in a region of densities between 0.01 and 0.1 nucleons fm -3 . A 
parameterization of the density dependent meson-nucleon coupling constants has been developed, which is based 
on microscopic Dirac Brueckner Hartree Fock (DBHF) calculations and reproduces the saturation point of nuclear 
matter as well as the bulk properties of finite nuclei. Since this parameterization yields a good description of normal 
nuclei, i.e. baryonic matter at normal density and small proton - neutron asymmetries but is also based on a realistic 
interaction which describes the scattering of two nucleons in the vacuum, it should give more reliable results for the 
baryonic structure at small densities and large proton-neutron asymmetries than corresponding calculations, which 
are based on purely phenomenological Hamiltonians like the Skyrme forces. The parameterization of the coupling 
constants in terms of analytic functions makes it easily accessible. 

The DDRMF calculations are performed for zero temperature as well as finite temperature in a periodic lattice of 
Wigner-Seitz (WS) cells of cubic shape. Pairing correlations are taken into account within the framework of the finite 
temperature BCS approximation assuming a density-dependent pairing force of zero range. 

At densities below 0.08 nucleons fm -3 the DDRMF approach predicts non-homogeneous structures, which are 
similar to those predicted by Thomas-Fermi or Skyrme Hartree-Fock calculations for the "pasta-phase" of matter in 
the crust of neutron stars. The structures obtained in the DDRMF calculations, however, are typically less pronounced 
than the corresponding structures resulting from Skyrme Hartree-Fock calculations. This may be a consequence of 
the finite range of the interaction in DDRMF as compared to the zero-range Skyrme forces. Also it turns out that 
the gain in energy due to the formation of non-homogeneous density distributions is slightly smaller in the DDRMF 
as compared to Skyrme HF. Also the DDRMF yields a smaller symmetry energy at the low densities, which provides 
smaller proton abundances. 

The effects of finite temperature tend to reduce the disposition of the baryonic matter to form non-homogeneous 
structures. This is driven by the larger entropy of the non-homogeneous matter as compared to the homogeneous 
phase. This leads to a reduction of the critical density for the formation of non-homogeneous structures. It drops 
from 0.08 nucleons fm~ 3 at T — to 0.05 nucleons fm -3 at T = 5 MeV. At densities below the critical densities, 
the gain in the free energy of the non-homogeneous as compared to the homogeneous realization is smaller and the 
resulting geometrical structures are less pronounced in the case of finite temperature. 

An attempt has been made to reproduce the results of the microscopic Skyrme HF and DDRMF calculation 
by means of the Thomas-Fermi (TF) approximation. Main features like the critical densities for the formation of 
geometrical structures of various shapes, the energies and the proton abundances can be reproduced, if one considers 
a parameter for the surface term, which is reduced with increasing proton asymmetry and with temperature. 

This work has been supported by the European Graduate School "Hadrons in Vacuum in Nuclei and Stars" (Basel, 
Graz, Tubingen), which obtains financial support by the DFG and received support from the GSI project TUEMUE 
and by the Spanish Ministry of Education and Science under grant no. SB-2005-0131. 
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